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Abstract 

We compare the dynamics of a Bose-Einstein condensate in two coupled 
potential wells with atoms diffracting from a standing light wave. The cor- 
responding Hamiltonians have an identical appearance, but with a different 
set of commutation rules. Well-known diffraction phenomena as Pendellosung 
oscillations between opposite momenta in the case of Bragg diffraction, and 
adiabatic transitions between momentum states are shown to have analogies 
in the two-well case. They represent the collective exchange of a fixed number 
of atoms between the wells. 
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I. INTRODUCTION 



The most common approach to the description of a trapped Bose-condensed gas is based 
on the mean-field approximation, which yields the Gross-Pitaevski equation for the macro- 
scopic wave function. This wave function, which depends on the number of atoms, plays the 
role of the mode function for the Maxwell field. This approach is reliable when the conden- 
sate is trapped in a single quantum state in a potential well. However, when the condensate 
is separated into two or more parts, so that more than one quantum state is populated, the 
mean-field approach is not evidently justified. It has been shown by Javanainen and Yoo 
that two originally separate parts of a condensate that are initially in a Fock state and that 
are brought to overlap will reveal an interference pattern that varies in position from one 
realization to another. This effect, which has also been observed experimentally [0] cannot 
be described by a single macroscopic wavefunction. A simple model for a condensate in a 
double potential well is defined by a field-theoretical Hamiltonian for a boson-Hubbard dimer 
0,^, which can be expressed in terms of SU(2) angular momentum-type operators with a 
quadratic term. This latter term represents the interaction between atoms in a well. The 
mean-field approximation is basically equivalent to classical equations of motion for the ex- 
pectation values of the SU(2) operators 0,3. The quantum regime has mainly been studied 
numerically, leading to collapse and revival [^], and to non-classical dynamics arising from 
the periodic modulation of the coupling between the wells 0. The formation of a two- well 
condensate by the raising of the barrier has been analyzed theoretically [§]. The situation of 
a Bose-Einstein condensate (BEG) in a two-well trap is also studied experimentally [§,0. 

A very similar Hamiltonian describes the situation of an atom diffracting from a standing- 
wave optical potential. This problem has received attention already in the early days of laser 
cooling More recent work has developed the band structure of the energy spectrum 



12 1, and a number of regimes have been distinguished that allow an analytical description 



p!3| . In a simple version of the model, the Hamiltonian is identical in form as in the two- well 



problem mentioned above. Now the quadratic term represent the kinetic energy of the atom. 
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The only difference between the two cases is that the commutation rules for the operators 
in the diffraction case are slightly simplified compared to the case of SU(2) symmetry. 

In this paper we discuss the analogy and the differences between these two systems. We 
point that a number of analytical solutions known for the diffraction problem can be carried 
over to the two-well system. The physics of these cases is discussed. 



II. BEC IN A DOUBLE POTENTIAL WELL 

We consider a potential consisting of two wells. When the barrier between the wells is 
not too low, the ground state and the first excited state \g) and |e) of a single atom are 
well approximated as the even and odd superposition of the lowest bound states in the two 
wells. Therefore, these states can be described as 

b)-^(|i) + |2»; |e) = i=(|l)-|2)), (1) 

with |1) and |2) the locahzed states in either well. When the energy separation between the 
excited and the ground state is indicated as h5, the off-diagonal element of the one-particle 
Hamiltonian Hi between the localized states is 

|2) = -h6/2 . 

At the low energies that are of interest here, the two particle interaction is well approximated 
by the standard contact potential U (r, r') — {4:nh^a/m)S{f—'r'), with a the scattering length. 
The second-quantized field operator is now 

in terms of the wavefunctions ■^j and the annihilation operators Sj of the single-particle states. 
The annihilation operators and the corresponding creation operators obey the standard 
bosonic commutation rules. The corresponding Hamiltonian is 



The wavefunctions ipi and 1^2 of the locahzed states have the same form, and we assume 
that they do not overlap. Then the interaction term can be expressed exclusively in the 
parameter k defined by 



hK= dr\i)i{r)W (4) 



m 

which measures the strength of the interatomic interaction. Performing the integrations in 
eq. d^) leads to the expression for the Hamiltonian 

-f/" = (aia2 + aafli) + ^ (aifliaiai + 02222222) , (5) 

where we took the zero of energy halfway the two energy levels of a single atom. This is 
also known as the boson- Hubbard dimer Hamiltonian ||^. 

The Hamiltonian (j^) can also be expressed in terms of SU(2) operators by applying the 
standard Schwinger representation of two modes. This leads to the definition 

Jq = 2 (^iSi — ^222) ; J+ = 0)10,2 , J- = a^fli ■ (6) 

These operators are related to the Cartesian components of angular momentum by the 
standard relations J± = ^ iJy, and Jq = J^. They obey the commutation rules for 
angular momentum operators 

[Jo,J±] = ±J±, [J+,J-] = 2Jo, (7) 
which generate the SU(2) algebra. The Hamiltonian (^) can be rewritten in the form 



H = - y(J+ + J-) + riK^ + ^[N'- 2N) , (8) 

with = a\ai + 0,20.2 the operator for the total number of particles. Obviously, the Hamil- 
tonian commutes with A^, and it is block-diagonal in the number of particles A^. For 
each value of N, the Hamiltonian can be expressed as 



with the A^-particle Hamiltonian 



(9) 



where the operators are now restricted to the + 1 Fock states \n,N ~ n) with n = 
0,1,. . .N, with n particles in well 1, and N — n particles in well 2. In the language of 
angular momentum, this manifold of states corresponds to the angular-momentum quan- 
tum number J = N/2, and the 2 J + 1 Fock states are eigenstates of Jq with eigenvalue 
fi = n — N/2, with fi = —J,—J-\-l,...,J. Note that fi is half the difference of the particle 
number in the two wells. For an even number of particles, the angular-momentum quantum 
number J as well as the 'magnetic' quantum numbers are integer, whereas these number are 
half-integer in case of an odd number of particles. The action of the operators Jq and J± on 
the Fock states has the well-known behavior 



with f^ = y ( J + /i) ( J — /i + 1). The /^-dependence of the strength of the hopping operators 
J± reflects the bosonic accumulation factor, which favors the arrival of an additional bosonic 
atom in an already occupied state. 

When the quadratic term in eq. (P) would be replaced by a linear term, the evolution 
would be a uniform rotation in the 2 J + 1-dimensional state space with angular frequency 
v5M-/?. The presence of the quadratic term makes the dynamics considerably more com- 
plex. Therefore we compare this dynamics with another well-known case in which a similar 
quadratic term appears. 

III. STANDING- WAVE DIFFRACTION OF ATOMS 

The translational motion of a two-level atom in a far detuned standing-wave light field 
is described by the effective Hamiltonian 



Jo l/i) = l/i) , J+ 1/^) = l/i + 1) , J- l/i) = 1/^ - 1) 



(10) 





A 




(11) 
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with ^ = ojo — oj the difference of the resonance frequency and the optical frequency, and 
ujr the Rabi frequency of each of the travelhng waves that make up the standing wave. The 
Hamiltonian takes a particularly simple form in momentum representation, since the kinetic- 
energy term is diagonal in momentum, and the potential energy changes the momentum 
by ±2hk. Therefore, we introduce momentum eigenstates which have the momentum 
2^hk. Then apart from an irrelevant constant, the Hamiltonian ([Tl|) can be represented in 
the algebraic form 

Hd = -— (5+ + 5_) + hKBl , (12) 

where k = 2hk'^/m determines the kinetic energy term, and 5 = u;|j/2A the atom-field 
coupling. The operators occurring on the r.h.s. are defined by the relations 

Bo\fi) = fi\fi); B±\fi) = \fi±l) . (13) 

They differ from the corresponding relations (|T0|) in that now the strength of the hopping 
operators is uniform. 

This Hamiltonian (0) has the same form as eq. (^, even though they describe com- 
pletely different physical situations. The difference is mathematically characterized by the 
commutation relations. The SU(2) relations (^ are replaced by the simpler set 

[5o,fi±] = ±fi±, [B+,B_]=0, (14) 

which is easily found from their explicit expressions (p!3D. The two operators B± are found 
to commute. A result of this difference is that the state space in the two-well case has a 
finite dimension 2 J + 1 = + 1, whereas the momentum space has an infinite number of 
dimensions. 

A mathematically identical set of operators occurs in the description of the dynamics of 
the Wannier-Stark system, consisting of a particle in a periodic potential with an additional 



uniform force WM. In that case, the eigenstates of Bq represent the spatially localized 



Wannier states, rather than the momentum states. 



We recall three approximate solutions of the evolution governed by the Hamiltonian (f2|), 
which are valid in different situations, and which allow analytical solutions. 

The Raman-Nath regime is valid for interaction times that are so short that the atom 



has no time to propagate. Then the quadratic term in (|T^) can be neglected, and the 
evolution is determined by the atom-field coupling S{t). The evolution operator is simply 
U = exp[i0(i?+ + i?_)/2], where = / dt6(t) is the integral of the coupling constant over 
the evolution period. The matrix elements of the resulting evolution operator for the pulse 



can be found by operator algebra in the form |l4 



{^,'\u\^l) = ^^'-^J,,.M (15) 

in terms of Bessel functions. For an initial state with a well-determined momentum, the 
time-dependent state following the pulse can be expressed as 

|vl/(t))^5:e-*'^"|/.')(/.'|f/|/i). (16) 



This leads to explicit analytical expressions for diffraction experiments [O]. The probability 
of transfer of n units of momentum is proportional to | J„(0)p. 

The Bragg regime is valid when the coupling 6 between neighboring momentum states 
is small compared to the kinetic-energy separation ^ of the initial state from its 

neighboring states |yU + 1) . This initial state leads to an oscillating time- dependent state 
between the two states and | — /i) with the same kinetic energy 

|vl/(t)) = cos^ |/x)+^sin^ h/i) , (17) 

apart from an overall phase factor. This can only occur when the momentum transfer 2/^ 
(in units of 2hk) is an integer, which corresponds precisely to the Bragg condition. The 
Pendelldsung frequency is given by fi^ = 5{5 /2kY^^^ / [(2/i — 1)!]^ ||13[. This expression is 
fully analogous to the effective Rabi frequency for a resonant multiphoton transition, with 
non- resonant intermediate states |T^JT^. 
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The regime of adiabatic coupling arises for a time-dependent atom-field coupling S{t) 
that varies sufficiently slowly, so that an initial energy eigenstate remains an eigenstate. 
The adiabaticity condition in the present case reads 

- « (18) 

When an atom passes a standing wave with a sufficiently smooth variation of the intensity, 
and the Bragg condition is fulfilled, the presence of two initially degenerate eigenstates | ± /x) 
leads to interference after the passage, which produces two outgoing beams. Because of the 
similarity between the two Hamiltonians (P) and ([T^), these well-known diffraction cases can 
be expected to have analogies in the dynamics of the two-well problem. 



IV. SYMMETRY CONSIDERATIONS OF GENERIC HAMILTONIAN 

The Hamiltonians and (0) can be represented in the generic form 

H = -h5L^ + HkLI, (19) 

with Lx = {L^ + L_)/2, = Lq, where the operators Lj represent Jj or Bi, depending 
on the commutation rules and the corresponding algebra that they obey. In the two-well 
case, the eigenstates \fi) of the operator represent number states in the two-well case, 
with the eigenvalue fi half the number difference between the wells. In the diffraction 
case, the states \fi) are momentum eigenstates. In this latter case, the coupling between 
neighboring momentum states is independent of fi (eq. (0)), whereas in the two- well case the 
/i-dependence of the hopping operator indicated in eq. ([T0|) reflects the bosonic accumulation 
effect. A consequence of this is also that the Hamiltonian in the diffraction case couples an 
infinite number of states \fi), whereas in the two-well case the number of coupled states has 
the finite value + 1. In the diffraction case we restrict ourselves to the situation that the 
Bragg condition is respected. Therefore, both in the diffraction case and in the two-well case 
II attains either integer or half-integer values. The action of Lz is the same in both cases. 
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The Hamiltonian (|T9D is invariant for inversion of /i. In order to demonstrate this, 
we introduce the inversion operator P, defined by the relation = | — //). In the 

diffraction case, the operator P corresponds to inversion of momentum, which does not 
change the kinetic energy. In the two- well case, the operator P represents interchanging 
the particle numbers in the two wells, which has no effect on the interparticle interaction. 
The commutation rules of the inversion operator with the operators Li are specified by 
PLzP = —Lz, PL±P = L=p, so that P inverts Ly and Lz, and commutes with L^. It 



follows that the Hamiltonian (|l^) commutes with P, so that it is invariant for inversion 
of /i. Therefore the Hamiltonian has vanishing matrix elements between the even and the 
odd subspaces, which are the eigenspaces of P with eigenvalue 1 and —1 respectively. For 
half-integer /x-values, these spaces are spanned by the states 

1/^)+ = — ^ — ; = — ^ — ; (20) 

for positive values of /i. In the case of integer yU- values, the state |/i = 0) also belongs to the 
even subspace. The even and odd subspace evolve independently from one another. This 
symmetry property of H depends on the fact that it is quadratic in the operator Lz- 

The action of the quadratic term in the Hamiltonian ( ]TP| ) on the new basis is simply 
given by the relation = /i^|/i)±. The action of the coupling term in the Hamiltonian 

can be expressed in a general form by introducing coefficients P^ for non-negative values 
of /i. In the case of the SU(2) algebra, we define P^ = whereas in the diffraction case 
we simply have P^ = 1. The matrix elements of can be fully expressed in terms of the 
coefficients P^ for positive /i. Within the even or the odd subspace, the operator has 
off-diagonal matrix elements only between two states for which the values of /i differ by 1, 
and we find 



12)^ = -F^^,, (21) 

provided that the value of fi is positive. These matrix elements coincide with those on 
the basis of the states For the state |/i = 0), which belongs to the even subspace of a 
manifold of states with integer /x-values, the matrix element is 



+ 



(i|l,|o) = Fi/v^. 



(22) 



On the other hand, in a manifold of states with half-integer /x- values, has a single non-zero 



diagonal element for /x = 1/2, that is given by 



±(l/2|L,|l/2)± = ±Fi/2. 



(23) 



Hence, in the case of half-integer //-values, the Hamiltonian projected on the even and the 



For integer values of the Hamiltonian for the odd subspace is identical to the Hamiltonian 

for the even subspace with fi )p 1. The only difference is that the even subspace also contains 
the state |0), which is coupled to the other states by the matrix element 



In both cases, the difference between the Hamiltonian parts on the even and odd sub- 
spaces are proportional to S. These differences are responsible for the energy splitting 
between the even and the odd energy eigenstates. Moreover, since these differences in the 
Hamiltonian parts occur for low values of /i, we expect that the even-odd energy splittings 
get small for large //-values. This is confirmed by numerical calculations. In Figs. 1 and 
2 we display the energy levels of the Hamiltonian, for a few values of 5/k, both for the 
double- well case (with = 100), and for the diffraction case. The energy levels are found 
to be alternatingly even and odd, with increasing energy. In the two-well case, the energy 
shifts and splittings due to the coupling are larger for the same value oi 5/k and the same 
value of /i. This arises from the factor F^, which is unity in the diffraction case, whereas 
in the two- well case it decreases from ~ J = A^/2 at /t = to zero at /t = J. In fact, the 
condition for weak coupling is that matrix elements coupling the states |/t) and \ii — 1) are 
small compared with their unperturbed energy separation. This condition can be expressed 
as 



odd subspace differ exclusively in the diagonal matrix element for /i = for which we find 





(25) 



10 



In the two- well case, the lowest energy states start out to be nearly equidistant up to /x-values 
where approaches one. 

V. PENDELLOSUNG OSCILLATIONS 

The energy splittings between the even and the odd eigenstates give rise to time- 
dependent states that oscillate between the states |±/i). In the diffraction case, they corre- 
spond to the well-known Pendellosung oscillations in the Bragg regime. Here we show that 
similar oscillations can occur for the two-well problem, and we give an analytical estimation 
of the oscillation frequencies. For the generic Hamiltonian given by (|l^), the Bragg condition 



is fulfilled when the inequality (|26|) holds. 

The energy differences between the even and odd states to lowest order in can be found 

from the effective Hamiltonian for two degenerate states that are coupled via a number of 
non-resonant intermediate states. This situation occurs for the states |±yu), with their 2yU — 1 
intermediate states. In this case, the intermediate states can be eliminated adiabatically, as 



demonstrated in Sec. 18.7 of ref. The resulting effective Hamiltonian for these two states 



|±/i) has an off-diagonal element that is the ratio between two products. The numerator 
contains the product of the successive 2/i matrix elements —h5F^j_i/2 of the Hamiltonian 
coupling neighboring states, and the denominator is the product of the 2/i — 1 unperturbed 
energy differences of the degenerate states |±/i) with the successive intermediate states. In 



the diffraction case, this result coincides with the calculation given in ref. [|T^], which was 
obtained by diagonalizing a tridiagonal matrix and keeping only the lowest order in 5/ n. 

Generalizing this result to the present case of the two states |±/i), we find that the 
effective Hamiltonian has the diagonal element 

{±li\%ff\±li) = hKii'' (27) 

and the off-diagonal element 
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{Tfi\H,jf\±fi) = -hn,/2, (28) 
with Qfj^ an effective oscillation frequency given by 

2+1 1 1 

The factor F is just the product of the coefficients successively coupling the states 
intermediate between \fi) and |— yu). In the diffraction case, we simply have F = 1, whereas 
in the case of SU(2) symmetry, applying to the double well we find 

These expressions are valid both for integer and half-integer values of /i. The eigenstates 
of the effective Hamiltonian are the even and odd states, and the eigenvalue equations are 
Heff = {fiKu'^ =1= For integer values of /i, the frequency fi^ is negative, so 

that the even states are shifted upwards, and the odd states are shifted downwards in 
energy. The opposite is true for half-integer values of ^. In both cases, the ground state is 
even, and the energy eigenstates for increasing energy are alternatingly even and odd. In 
view of the results of the numerical calculation mentioned above, one may expect that this 
alternating behavior of the even and odd eigenstates is valid for all finite values of the ratio 
6/k. 

For an initial state this effective Hamiltonian leads to a time-dependent state that is 
given by ([TTD , apart from an irrelevant overall phase factor. This shows that the oscillating 
solution ([I7|) corresponding to the Bragg regime of diffraction can be generalized to the case 
of a condensate in a double well. The same expression (|1^ remains valid, while the oscillation 
frequency fl^ is determined by eqs. p9| ) and (pOD. This describes a state of the condensate 
atoms in the double well in the weak-coupling limit. In this case, the state oscillates between 
the Fock states ria) = \N/2 + /i, N/2 - fi) and ns) = \N/2 - /i, N/2 + fi). 

Obviously, when the initial state is given by the system is in a stationary state, 

and no oscillations occur. In this case, Pendellosung oscillations can still be induced by 
including in the Hamiltonian a term that is linear in L^. In the diffraction case, there is 
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no obvious physical realization of such a term. For the Wannier-Stark system, where the 
quadratic term in is absent, the linear term can be realized by imposing a uniform force, 



which gives rise to Bloch oscillations |14| , pT| . In the case of the BEC in a double well a term 
hE,Lz/2 in the Hamiltonian can be realized by imposing an energy difference TiS, between the 
single-particle ground states in the two wells. When this term is periodically varying, it can 
be used for coherent control of the condensate The additional term couples the even 



and odd subspaces, thereby breaking the symmetry opf the Hamiltonian. On the basis of 
the states the effective Hamiltonian attains the off-diagonal element 

^{y^\%ff\^)^ = K^^/2. (31) 

When we assume that both 5 and ^ are small compared with the splitting due to the 
interparticle interaction k, so that we remain in the Bragg regime, the two states |±/i) 
remain decoupled from the other number states, and we have an effective two-state system. 
In practice, the parameter can be easily controlled, so that many effects of two-state 
atoms [113 can also realized for these two states. For example, in analogy to the excitation 



of ground-state by an adiabatic sweep across the resonance, one could create an effective 
transfer from the state to the state |— /u) by varying the parameter ^ adiabatically from 
a positive to a negative value that is large compared to Vt^. This gives an effective collective 
transfer of n = 2fi atoms from one well to the other one. 

VI. TIME-DEPENDENT COUPLING 

When the coupling 6{t) varies with time, the time-dependent eigenstates of the Hamil- 
tonian are coupled to each other. The eigenstate that correlates in the limit 5 to the 
state l/i) , is denoted as ^t)- Note that even eigenstates are only coupled to other even 



eigenstates, and odd eigenstates to odd eigenstates. The coupling results from the time 
dependence of the eigenstates. In fact, the term in the Schrodinger equation coupling 
to \ip^) is proportional to 
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it) I = - {v^ (t)| 4 I'fit it)) E±-l± ^ ^^"^ (^2) 



This couphng is ineffective in the case that the r.h.s. of eq. (^) is small compared with 
[E^^ — E'^)/h. In this case, an initial eigenstate remains an eigenstate at all times. This is 



the standard case of adiabatic following, which has been discussed in the diffraction case |]T3 
Since within the even or the odd subspace there are no degeneracies, the dynamics of adia- 
batic following is particularly simple. When the coupling coefficient 6 is smoothly switched 
on, with the system initially in the state = (1/^)^ + |/i)_)/v^, the time-dependent state 
is obviously 



with 

i^{t) = J'dt' (E+{t') + E^it')) /2h the average phase, and r]{t) = j' dt' [E+{t') - E^{t')) /U 
the accumulated phase difference of the two eigenstates. In a time interval that the coupling 
6 is constant, the phase difference r]{t) increases linearly with time, and the state (|33| ) gives 
rise to expectation values oscillating at the single frequency (^E^{t') — E^{t')^ /h. When 
the coupling is switched off again, the phase difference approaches a constant limiting value 
rj = T] (oo). The state (P^D at later times corresponds to a linear superposition of the states 
|±/i) proportional to cos(r//2) + \ sin(r//2). Again, as we see, this effect that is known 
in the diffraction case also has a counterpart for the double-well problem, where adiabatic 
switching of the coupling between the wells leads to a linear superposition of the Fock states 
["■i, ^2) = \N/2 + fi, N/2 — jj) and |ni, 712) = \N/2 — /i, N/2 + /i). By proper tailoring of the 
pulse, the final state can be made to coincide with either one of these Fock states, or with 
the even state or with the odd state 

In contrast, when the coupling term 5{t) has the form of a short pulse around time 0, 
such that the action of the quadratic term can be neglected during the pulse, the initial 
state l/i) couples to all other states . The state vector has exactly the same form (|16[) as 
for diffraction in the Raman-Nath regime. For the two-well problem, the evolution operator 
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takes the form U = exp{i(l)Lx), with (f) = J dt5{t), which has matrix elements that can be 
expressed in the Wigner rotation matrices by 

(/i'|[/|/i)=z^'-^4^(0), (34) 

with J = N/2. A comparison with Eq. (|1^) shows that for the two- well-problem, the Wigner 
functions play the same role as the Bessel functions in the diffraction case. 

In Fig. 3 we plot the energy difference E'^ — E~ between the even and odd eigenstate 
in the two-well case, for = 100, and for a few values of 6/k,. The splittings decrease 
monotonously for increasing quantum number /i. 



VII. CONCLUSION 

In this paper we have analyzed both the similarity and the difference between the dy- 
namical behavior of atom diffraction from a standing wave and a Bose-Einstein condensate 
in a double- well potential. In both cases, the Hamiltonian is given by the generic form (|19|), 



the only difference being in the commutation rules for the operators Lj. with i = x,y, z. 
Well-known diffraction phenomena as Pendellosung oscillations between opposite momenta 
in the case of Bragg diffraction, and the result of adiabatic transitions between momentum 
states have counterparts in the behavior of the atom distribution over the two wells, in the 
case that the coupling between the wells is weak compared to the interatomic interaction 
or slowly varying with time. A common underlying reason for these effects is the symmetry 
of the Hamiltonian for inversion fi ^ —fi, and the energy splitting between even and odd 
states arising from the coupling term. In these cases, effective coupling occurs between the 
states \ni,n2) and 1^2, ni) with opposite imbalance between the particle numbers in the two 
wells. These states are coupled without population of the intermediate states, so that a 
number of rii — n2 particles oscillate collectively between the two wells. The interparticle 
interaction is essential for this effect to occur. A simple analytical expression is obtained 
for the Pendellosung frequency. An initial state |ni,n2) with a well-determined number of 
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atoms in each well can be transferred to a linear superposition of |ni, ^2) and \n2, rii), which 
is a highly entangled state of the two wells. A similar analogy is obtained to diffraction 
in the Raman-Nath regime. For the double-well problem this requires that the coupling is 
sufficiently short to ignore dynamical effect of the atomic interaction during the coupling. 
The well-known diffraction pattern in terms of the Bessel function is replaced by elements 
of theWigner rotation matrix for the double well. These effects do not show up in the 
mean-field approximation, where the Gross-Pitaevski equation holds. 
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FIG. 1. Energy levels in units of Hk for the double well with = 100 particles, for various 
values of — . The levels are labeled by the quantum number /i. 
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FIG. 2. Energy levels in units of fin for the diffraction case, for various values of 
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